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p ,! Abstract 
D 

l/^ ' The use of the hyperspherical harmonic (HH) basis in the description of bound states in an A- 

' 

body system composed by identical particles is normally preceded by a symmetrization procedure 
in which the statistic of the system is taken into account. This preliminary step is not strictly 
necessary; the direct use of the HH basis is possible, even if the basis has not a well defined behavior 

T3 ■ 

^ ' under particle permutations. In fact, after the diagonalization of the Hamiltonian matrix, the 

a : 

eigenvectors reflect the symmetries present in it. They have well defined symmetry under particle 
^ , permutation and the identification of the physical states is possible, as it will be shown in specific 

' cases. The problem related to the large degeneration of the basis is circumvented by constructing 

■ the Hamiltonian matrix as a sum of products of sparse matrices. This particular representation 

, of the Hamiltonian is well suited for a numerical iterative diagonalization, where only the action 

' of the matrix on a vector is needed. As an example we compute bound states for systems with 

^ = 3 — 6 particles interacting through a short-range central interaction. We also consider the 
case in which the potential is restricted to act in relative s-waves with and without the inclusion 
of the Coulomb potential. This very simple model predicts results in qualitative good agreement 
with the experimental data and it represents a first step in a project dedicated to the use of the 
HH basis to describe bound and low energy scattering states in light nuclei. 



PACS numbers: 31.15.xj, 03.65. Ge, 36.40.-c, 21.45.-v 



1 



I. INTRODUCTION 



The ab initio description of light nuclear systems, starting from the nucleon-nucleon (NN) 
interaction, requires well established methods to solve the Schrodinger equation. Among 
them, the Green function Monte Carlo (GFMC) method has been extensively used to de- 
scribe light nuclei up to A = 10 and the no-core shell model (NCSM) up to A = 12 l|, l2|. In 
the A < 4 systems, well established methods for treating both bound and scattering states 
exist as the Faddeev equations {A = 3) and the Faddeev-Yakubovsky equations {A = 4) in 
configuration or momentum space, and the Hyperspherical Harmonic (HH) expansion. All 
these methods have proven to be of great accuracy and they have been tested using different 
benchmarks 3|-|5|. 

The HH method provides a systematic way of constructing a complete basis for the 
expansion of the A-particle wave function and its use in the A > 4 systems has been subject 
of intense investigations over the last years. In the specific case of application to nuclear 
physics, the wave function has to be antisymmetric and, therefore, the HH basis has been 
managed to produce basis states having well defined properties under particle permutations. 
Different schemes to construct hyperspherical functions with an arbitrary permutational 



symmetry are given in Refs. 



m 



Recently, a procedure for constructing HH functions in 



terms of a single particle basis has been proposed in Ref. j^. 

In a different approach, the authors have used the HH basis, without a previous sym 



metrization procedure, to describe bound states in three- and four-particle systems |10 |. 
It has been observed that the eigenvectors of the Hamiltonian matrix reflects the symme- 
tries present in it, even if it has been constructed using the non-symmetrized basis. The 
only requirement was to include all the HH basis elements having the same grand angu- 
lar quantum number K. It is a property of the HH basis that basis elements having well 
defined behavior under particle permutation can be constructed as a linear combination of 
HH elements having the same value of K. Therefore, if the Hamiltonian commutes with the 
group of permutations of A objects, Sa, the diagonalization procedure generates eigenvectors 
having well defined permutation symmetry that can be organized in accordance with the 
irreducible representations of 5*^. Moreover, identifying those eigenvectors with the desired 
symmetry, the corresponding energies can be considered variational estimates. In particular. 



in Ref. 



10| , it was possible to identify a subset of eigenvectors and eigenvalues corresponding 



exactly to those that would be obtained performing the preliminary symmetrization of the 
states. It should be noticed that the simplicity of using the HH basis without a preliminary 
antisymmetrization step, has to be counterbalanced with the large dimension of the matrices 
to be diagonalized. However, at present, different techniques are available to treat (at least 
partially) this problem. 

In the present article we continue the study of the non-symmetrized HH basis, extending 
the applications to systems with A > 4. In pursuit of this goal, we have developed a 
particular representation of the Hamiltonian matrix, which is systematic with respect to 
the number of particles and well suited for a numerical implementation. As mentioned, one 
of the main problem in using the HH basis is its large degeneracy, resulting in very large 
matrices. On the other hand, the potential energy matrix, expressed as a sum of pairwise 
interactions, cannot connect arbitrary basis elements differing in some specific quantum 
numbers. This means that in some representation each pairwise-interaction term has to be 
represented by a sparse matrix. For example, the matrix representation of the potential 
V{1,2), constructed in terms of basis elements in which the quantum numbers of particles 
(1, 2) are well defined, is sparse in A > 3 systems. In fact, its matrix elements connecting 
basis elements with different quantum numbers labelling states which do not involve particles 
(1,2) are zero. A problem arises when the matrix elements of the generic term V{i,j), 
defining the interaction between particles has to be calculated using basis elements in 

which the quantum numbers of particles are not well defined. One operative way to 
solve this problem consists in rotating the basis to a system of coordinates in which particles 
have well defined quantum numbers. This makes the matrix V{i,j) sparse. However, 
we would like the rotation matrix to be sparse too, which in general it is not true. This 
last problem is solved noticing that the rotation matrix can be expressed as a product of 
sparse matrices, each one representing a rotation which involves a permutation of particles 
of successive numbering. After these manipulations the potential energy matrix results in a 
sum of products of sparse matrices suitable for numerical implementations. 

An advantage in using the non-symmetrized HH basis appears when symmetry breaking 
terms are present in the Hamiltonian. In the case of the nuclear Hamiltonian with charge- 
symmetry breaking terms, this means that different total isospin components T are present 
in the wave function. For example, the three-nucleon bound state wave function includes T = 
1/2,3/2 components and the four-nucleon bound state wave function includes T = 0, 1,2 
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components, requiring the inclusion of different spatial symmetries in the wave function. 
Therefore, considering all the possible spin and isospin components, the number of HH 
states having well defined spatial symmetries, necessary to construct the wave function, 
and the dimension of the non-symmetrized basis is comparable. High isospin components 
are in general a small part of the total wave function. They are difficult to include in the 
antisymmetrized basis since appreciably increases the number of basis elements and, at the 
same time, they improve very little the description of the state. In practical cases they are 
disregarded, or partially included, with the consequence that the occupation probabilities 



of the high isospin states are not always well determined (see Ref. Conversely, using 

the non-symmetrized basis, all the isospin components are automatically generated. As an 
example we will show results for A = 3 — 6 systems using short-range central interactions 
with and without the inclusion of the Coulomb potential. 

To summarize, in this paper we present the implementation of the non-symmetrized HH 
basis for A-hodj system using the factorization of the potential energy matrix mentioned 
before. In order to give a detailed description of this construction, we consider only spatial 
degrees of freedom; accordingly, we show examples using a central interaction. The diago- 
nalization of the Hamiltonian produces eigenvectors organized in multiplets of the dimension 
of the corresponding irreducible representation of Sa, and the different symmetries will be 
identified using the appropriate Casimir operator. Not all the states belonging to a partic- 
ular representation can be antisymmetrized using the spin-isospin functions of A nucleons 
and, therefore, these states are not physical. It should be noticed that the physical states 
could appear in very high positions of the spectrum, in particular this is the case for A > A 
systems. On the other hand, the iterative methods, as the Lanczos method, used to search 
selected eigenvalues and eigenvectors of large matrices are more efficient for the extreme 
ones. To this respect, we have found very convenient to use the symmetry-adapted Lanczos 



method proposed in Ref. 12[, which restricts the search to those states having a particular 
symmetry. When possible, comparisons to different results in the literature will be done. 
Since we have in mind the description of light nuclear systems using realistic interactions, 
this study can be considered a preliminary step in the use of this technique. 

The paper is organized as follows; section [III is devoted to a brief description of the HH 
basis. In sections IIIII the expression for the potential energy matrix in terms of HH states 
is given. In section [IV] the results for the models proposed are shown. Section |V] includes a 



brief discussion of the results and the perspectives of the present work. 



II. THE HARMONIC HYPERSPHERICAL BASIS FOR A BODIES 

In this section we introduce the notation and we present a brief overview of the properties 
of the HH basis. 

A. Basic properties of the HH basis 



In accord with Ref. 



10|, we start with the following definition of the Jacobi coordinates 



for an A body system with Cartesian coordinates ri . . . 



where m is a reference mass, N = A — 1, and where we have defined 

2=1 ^ i=l 

Let us note that if all the masses are equal, rrii = m , Eq. ([1]) simplifies to 



^N-j+i = ^1 (r^+i - Xj-) , J = 1, . . . , AT . (3) 
For a given set of Jacobi coordinates xi, . . . , x^r, we can introduce the hyperradius p 

/ ^ \ 1/2 / ^ \ 1/2 / 9 ^ \ 1/2 

/^=(E-n =(2E(r^-X)' =(lE(r.-r.)M , (4) 
and the hyperangular coordinates 

= {xi, . . . ,Xn,4>2, ■ ■ ■ ,4>n) , (5) 
with the hyperangles (pi defined via 

cos (f)i = , „ = , i = 2,...,N. (6) 

^^xi H \-xf 
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The radial components of the Jacobi coordinates can be expressed in terms of the hyper- 
spherical coordinates 

Xn = p cos 07V 

Xm-1 = psin</)7vCOS0Ar_i 



Xi = p sm ■ ■ ■ sm cos ( 



(7) 



X2 = p sm 0Ar ■ ■ ■ sm 03 COS 02 

Xi = p sin 0JV ■ ■ ■ sin 03 sin 02 . 
Using the above hyperspherical angles Qn, the surface element becomes 

TV 

dflN = sin 9i dOi d(pi JJ^ sin 6j dOj dip j {cos 0j)^(sin (j)j)^^~'^d(j)j , 
i=2 

and the Laplacian operator 

N 



\dp^ p dp P^ J 



(8) 



(9) 



where the A^(f27v) is the generalization of the angular momentum and is called grand angular 
operator. 

The HH functions 3^[a'](^7v) are the eigenvectors of the grand angular momentum operator 



A%{nr,) + K{K + 3N- 2) J 3^[x](fi,v) = . (10) 

They can be expressed in terms of the usual harmonic functions Yim{x) and of the Jacobi 
polynomials Pf^'^{z). In fact, the explicit expression for the HH functions is 



N) 



■ N 

.i=i 



■ N 

n '"^k 

J=2 



(11) 



where [K] stands for the set of quantum numbers li, . . . ■ ■ ■ ,^N,n2, ■ ■ ■ ,nN, and the 

hyperspherical polynomial is 

0)p%'"-.-i (^^.) ^ (cos0,)'^(sin 0,)^-^P.7-^'"'^ (cos 20,) , (12) 



with the Kj quantum numbers defined as 



(13) 



1=1 



The normalization factor is 



I 2i2n + a + (3 + l)n\T{n + a + (3 + 1) 
" V r(n + a + l)r(n + /3 + l) ' ^ ^ 

where, for the special choice of hyperangles given by Eq. ([7]), = Kj + 3j/2 — 1 and 
ai- = lj+1/2. The quantum number K = Kn is also known as the grand angular momentum. 
The HH functions are normalized 

j dQN(^lK']{^N)^ y[K]{^N) = SlKUK'] , (15) 

moreover, the HH basis is complete 

J2 (ymi^N)) yud^N) = s'^'-'in'j, - n^) . (i6) 

[K] ^ ^ 

With the above definitions, the HH functions do not have well defined total orbital angular 
momentum L and z-projection M. It is possible to construct HH functions having well 
defined values of LM by coupling the functions Yi.mj{xj). This can be achieved using 
different coupling schemes. Accordingly, we can define the following HH function 



■ N 

J=2 



^A^-l) |ljv_i ®^i]v (^A^) 



LM 

(17) 



having well defined values of LM, with the particular coupling scheme in which particles 
(1, 2) are coupled to L2, which in turns, with I3, is coupled to L3 and so on, generating N — 2 
intermediate Lj-values. The set of quantum numbers [K] includes the n2. ■ - Un indices of 
the Jacobi poljTiomials, the li . . .1^ angular momenta of the particles and the intermediate 
couplings L2 . . . Lat-i. 

In the definition of the hyperspherical coordinates in terms of the radial components of 



the Jacobi coordinates it is useful to introduce the hyperspherical tree structure [13|. For ex- 
ample, the particular choice of Eq. ([7]), in the coupling scheme of Eq. ( [T71) .corresponds to the 
one depicted in Fig. [H where we can also read the above-mentioned angular- momentum cou- 
pling scheme. However, other definitions are possible, and the corresponding hyperspherical 
functions can be related using the T-coefficients Q, Schematically, these coefficients 



relate the following tree structures 




Ni 



i4',_i,n,_i 



' ni-iiii-iKi 



ni_i=0 




(18) 



Ki, rii 

and play the same role of three-momenta recoupling as the 6j coefficients, but for the grand- 
angular momenta. Here Ki = Ki^i + + 2ni = i^Tj-i + 2ni The explicit definition of the 
coefficients is 



r 



rii-l 



rii 



2 



(19) 



2"''+°'-!+' r(n,_i + + l)r(ni_i + a^^, + 1) 



rfy (1 - y)°'.-i+"'-(i + y)"'.p„?--^"'» (y) pr;^-"'' (y) . 



) ^i — 1 



In this expression the value of the coefficients Uj depend on the value of the partial grand 
angular momentum Kj or partial angular momentum /j, which labels the node or the leave 
respectively, and on the topology of the hyperspherical tree. Having in mind that in a binary 
tree a node and its child nodes form a sub-binary tree with Ma^ nodes and leaves the 
coefficients read 



1 



a,- 



(20) 



Furthermore, the integral in Eq. (I19p can be rewritten as an hypergeometrical function 
using the following identity: 



T{a-r + n)T{P + n + l)T{p + m + l)r(r + 1) 



mini 



4-^3 



r(p + l)r(a - r)r(/3 + r + ^ + 2) 
-m, p + cr + m + l,r + l,r + l,r — a + 1 
p + l,(3 + T + n + 2,T-a-n+l 



(21) 



For the sake of completeness, we also report the notation we use for the recoupling of three 
angular momenta 



(22) 

where we have defined 

^L'"'t'iL- = (-l)''-'^''-'^''''^V2Li_i + 1 •^2L,_i + 1 J I. (23) 

\ h Li Li_i J 

Both T- and T-coefficients have particular relevance in the construction of HH functions 
with arbitrary permutational symmetry js], [igI . 



B. Rotation matrices between HH basis elements of different Jacobi coordinates 

The Jacobi coordinates explicitly depend on the way of numbering the A particles. In 
particular, for an equal mass system, we have selected a successive order in Eq. ([3]) starting 
from the definition of xa? = r2 — ri. In the following we will refer to this set as the reference 
Jacobi set. However, different choices are possibles, starting for example from x^r = — fj, 
with the related HH functions depending differently on the particle variables. In general, 
the Jacobi coordinates can be defined from a permutation {p = pi . . . p^} of the A particles, 

. . . Fp^, resulting in a re-definition of the Jacobi coordinates in which rj, on Eq. ([3]), is 
changed to rp.. The associated HH functions, 3^j^^(f2^), are still defined by Eq. (fT7|) . The 
explicit indication of the index p of the permutation allows to trace back the dependence 
on the particle variables. It is a general property of the HH basis that elements constructed 
using a permutation p in the arrangement of the particles can be expressed as a linear 
combination of HH basis elements defined using some other order, both having the same 
grand angular quantum number. In our case, we use the HH basis constructed with the 
reference Jacobi set to express bases constructed with other arrangements. Accordingly, the 
property reads 

= E c^mym i^N) , (24) 

IK'] 

where the sum runs over all quantum numbers compatible with the condition K = K'. 
As indicated, in the transformation the total angular momentum LM is conserved. For a 
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given number of particles, Nk denotes the number of HH functions having the same value of 
K. Consequently, the coefficients of the transformation C*[^jx'] form a matrix of dimension 



Nk X Nk- For A = ?> these matrix elements are the Raynal-Revai coefficients [17|, whose 
expression is explicitly known. For A > 3 the coefficients cannot be given in a close form, 
and a few methods have been derived for their calculations 0, [l8l20|. 



Here we are interested in a particular set of coefficients relating the reference HH basis 
to a basis in which the ordering of two adjacent particles have been transposed. It is easy 
to verify that there are A — 1 sets of Jacobi coordinates of this kind based on the following 
ordering of the particles: (ri, . . . , r^, r^-i), (ri, . . . , r^_i, r^), • • • , (ri, rs, ra, . . . , r^), 
(r2,ri, . . . ,r^). This last ordering results in a Jacobi set in which all the Jacobi vectors 
are equal to those of the reference set except the last one, xjy, which is now x'^ = ri — Y2- 
The other A — 2 orderings lead to — 1 Jacobi sets that differ, with respect to the original 
Jacobi set, in the definition of two Jacobi vectors. In fact, given the transposition between 
particles j, j + 1, only the Jacobi vectors Xj and Xj+i, with i = N — j + 1, are different. We 
label them x^ and x^_,_]^, and explicitly they are 

x'- \.VU±WEEilA^. , 

— : Xj -| ; Xj+1 

i L (25) 

, _ vu + iy-2{j + i) 1 

J J 

with i = 1, . . . , N — 1. The value i = 1 corresponds to the transposition of the pair (r^-i, r^), 
whereas the value i = N — 1 corresponds to the transposition of the pair (r2,r3). Let us 
call 3^[^(fijv) HH basis element constructed in terms of a set of Jacobi coordinates in 
which the i-th and i + 1-th Jacobi vectors are given from Eq.( [25!) with all the other vectors 
equal to the original ones (transposed basis). The case i = N corresponds to the special 
case, mentioned before, in which all the vectors are equal except x^v- The coefficients 

•^mK'] = I dn4y(ii{{n^)ry(i\{^r^) , (26) 

are the matrix elements of a matrix Af^'^ that allows to express the transposed HH basis 
elements in terms of the reference basis. They are a particular case of the general C^k]\k'] 
defined in Eq.( [24l) and, therefore, the total angular momentum as well as the grand angular 
quantum number K are conserved in the above integral {K = K'). The coefficients -A^^KjiK'] 



can be calculated analytical 



Revai matrix elements 



14 



using the T- and T- coupling coefficients and the Raynal- 



15| . In fact, we have seen that only two Jacobi coordinates 
10 



are changed in the construction of the transposed HH basis (see Eq. (l25i) ). If the two 
coordinates are directly coupled both in grand-angular and angular space, as is the case 
for the pair Xi and X2, corresponding to i = 1, the coefficient reduces to the Raynal-Revai 
coefficient. Explicitly, 

N 



A 



l,LM 
[K][K'] 



6 



KJ<' 



i=3 



(27) 



with 



LM 



ILM 



(28) 



whose analytic form has been given in Ref. 



17j. When 2 < z < - 1, we still have a 



transformation between only two Jacobi coordinates, and the coefficients read 



A'- 



i,LM 
[K][K'] 



'i-l i-1 



k=2 



(i)j^i-iKi-\,Li^\Kij^i 
li,l^,li+i,l',LiKi,L'.K'. 



N 



N 



.j=i+2 



k=i+l 



(29) 



where Lj^ = L and Kf^r = K. From the conservation of partial angular and grand angular 

xf + xf,i^ the matrices ^*l4 can be obtained from a 



momenta and the fact that xf + xf^-^ 
three-dimensional integral. However, as has been shown in Ref. 15|, they can be reduced to 
Raynal-Revai coefficients using the T- and T-coefficients to recouple the quantum numbers 
relative to the Jacobi variables Xj and Xj+i. The final expression is 

Li fi, 

(30) 

where Ki = li + /j+i + 2nj. Finally, the case i = N corresponds to the transposition of 
particles (1, 2) resulting in x'^ = — x^r and the coefficient reduces to a simple phase factor 



aN,LM 
■^\K][K'] 



[Kim ■ 



(31) 



We are now interested in obtaining the rotation coefficients between the reference HH 
basis and a basis in which the last Jacobi vector is defined as x'^ = Vj — r^, without loosing 
generality we consider j > i. A generic rotation coefficient of this kind can be constructed 
as successive products of the A'^k][k'] coefficients. For j > 3, if XAr_j_|_2 is the last Jacobi 
vector in which particle j appears, at maximum 2(j — 2) factors have to be included in the 
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product. To see this, we start discussing the case j = 3 resulting in two different bases, 
one having the vector x'^ = ra — ri and the other the vector x'^ = — r2. The rotation 
coefficient that relates the basis having the vector x^ = r3 — ri to the reference basis, is 
'^[K][K']^^ since it corresponds to the transposition of particles (2,3). In the second case we 
have to consider the transpositions of particles (2,3) and (1,2) and, the coefficients results 
the same as before times the phase given by the coefficient A^^^^x']- Therefore, in the case 
j = 3, the rotation coefficients includes at maximum the multiplication of two ^-coefficients. 
There are three vectors x^ = rj — with j = 4. When i = 1 or i = 2, the transposition 
(3, 4) leads to the previous case and two and three factors are needed respectively. For the 
case x'^ = r4 — r3, the intermediate transposition (2, 3) is needed with the consequence that 
the rotation coefficient includes four factors, and so on. 

Defining y^j^^{Q^^) the HH basis element constructed in terms of a set of Jacobi coor- 
dinates in which the A^-th Jacobi vector is defined x^ = rj — r^, the rotation coefficient 
relating this basis to the reference basis can be given in the following form 

= / dn[y(i,'{in^)Yy(i,'{{n%) = [a^^^^ ■ ■ ■ a^^^] j^^^^,, . (32) 

The particular values of the indices ii, . . . ,in, labelling the matrices Af^^, . . . , Af^ , depend 
on the pair (i, j). The number of factors cannot be greater than 2[j — 2) and it increases, 
at maximum, by two units from j to j ' + 1. The matrix 

Blf = A^...A^, (33) 

is written as a product of the sparse matrices ^f^'^'s, a property which is particularly well 
suited for a numerical implementation of the potential energy matrix as is discussed in the 
next section. 

III. THE POTENTIAL ENERGY MATRIX IN TERMS OF THE ^-COEFFICIENTS 

The potential energy of an A-body system constructed in terms of two-body interactions 
reads 

V = Y,yihj) ■ (34) 

i<j 

12 



Considering the case of a central two-body interaction, its matrix elements in terms of the 
HH basis of Eq. f|T71) are 

VmK'ip) = J2{y(k'({^N)\vit,j)\y(kf]{n^)) . (35) 

In each element {y^j^ \ V {i, j)\yl"j^^) the integral is understood on all the hyperangular vari- 
ables and depends parametrically on p. Explicitly, for the pair (1,2), it results 

VmUip) = {yiK'{{^N)\v{i,2)\y(}i^{nr,)) = 

Sh,i[ ■ ■ ■ Si^^i'^6l2,l'2 ■ ■ ■ ^Ln,l'^^K2,k'^ ■ ■ ■ Skn,k'^ (36) 

The above formula shows that for A > 2 the matrix representation of V(l,2) is sparse in 
this basis. Using the rotation coefficients, a general term of the potential V{i,j) results 

VmK'iiP^ = {y[K]\^N)\v{z,j)\y;},^{n^)) = 

[K"] [K'"] 

It should be noticed that 

(3^[iT(^^)i^(i'2)i3^f4(fi^)) = {ymm\vi^,J)\ytK'f]i^i,)) , m 

therefore Eq. (l37Il results 

VmK']ip) = {y[K{i^N)\v{z,j)\y(},^\{n^)) = 

(39) 

'^[K"][K]'^IK"'][K'] ^lK"][K"']yPj ' 

[!<"] [/<"'] 

or, in matrix notation, 

VM = [B^''YVu{p)Bl^''. (40) 
The complete potential matrix energy results 

J2v,,{p) = J2[B^''fVMB^'^. (41) 

ij ij 

The matrices B^j'^^ are block matrices with each block labelled by the grand angular mo- 
mentum K. Moreover, each block is constructed as a product of the sparse matrices Af^ 
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as defined in Eq. (l33|) . On the other hand the matrix Vi2{p), defined in Eq.( l36|) . couples 
different values of K but it is diagonal in the quantum numbers related to particles 3, . . . ,A. 

Each term of the sum in Eq. fl41l) results in a product of sparse matrices, a property which 
allows an efficient implementation of matrix- vector product, key ingredient in the solution 
of the Schrodinger equation using iterative methods. 



IV. RESULTS FOR A = 3 - 6 SYSTEMS 



In this section we present results for A = 3—6 systems obtained by a direct diagonalization 
of the Hamiltonian of the system. The corresponding Hamiltonian matrix is obtained using 
the following orthonormal basis 



{pn\m[K]) = 



mi 



(42) 



(a + m)! 

where Lm\(3p) is a Laguerre polynomial with a = 3N — 1 and /3 a variational non-linear 
parameter. The matrix elements of the Hamiltonian are obtained after integrations in the 
p,Q spaces. They depend on the indices m^m' and [i^], [K'] as follows 



(m' [K']\H\m[K]) 



(T^n'm - K{K + 3iV - 2)T^:=4) 



(2) 



E 



LM rnj,LM -irm.m' 
[K][K"]'^[K"'][K'] ^[K"][K"'] 



[K"][K'"] 



(43) 



The matrices T^^^ and T^^^ have an analytical form and are given in Ref. 10|. The matrix 
elements V^'-^^,-^ are obtained after integrating the matrix Vi2{p) in p-space (we will call the 
corresponding matrix Vu)- Introducing the diagonal matrix D such that {[K'] \ D \ [K]) = 
^[K],[K']K{K+3N — 2), and the identity matrix I in i^-space, we can rewrite the Hamiltonian 
schematically as 

(44) 



H = _:^((i)T ® I + (2)r ® D) + YiB^'^Y B. 

m ^ ■' 



LM 
ij ) 



in which the tensor product character of the kinetic energy is explicitly given. A scheme to 
diagonalize such a matrix is given in the Appendix. 
We choose as central potential the Volkov potential 



V{r) = Vfie-'^'/^i + lOie-'^'/^' 



(45) 
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with Vr = 144.86 MeV, Ri = 0.82 fm, Va = -83.34 MeV, and R2 = 1.6 fm. The nu- 
cleons are considered to have the same mass chosen to be equal to the reference mass m 
and corresponding to h?/m = 41.47 MeVfm"^. With this parametrization of the potential, 
the two-nucleon system has a binding energy E2N = 0.54592 MeV and a scattering length 



= 10.082 fm. This potential has been used severaj 
use very useful to compare different methods 



11 



21 



times in the literature making its 



221 ]. The use of central potentials 



in general produces too much binding, in particular the A = 5 system results bounded. 
Conversely, the use of the s-wave version of the potential produces a spectrum much closer 
to the experimental situation. This is a direct consequence of the weakness of the nuclear 
interaction in p-waves. Accordingly, we analyze both versions of the potential, the central 
Volkov potential and the s-wave projected potential. The results are obtained after a direct 
diagonalization of the Hamiltonian matrix of Eq.f H3|) including rrimax + 1 Laguerre polyno- 
mials with a fix value of /3, and all HH states corresponding to maximum value of the grand 
angular momentum Kmax- The scale parameter /3 can be used as a non-linear parameter 
to study the convergence in the index m = 0, 1, . . . , rrimax, with rrimax the maximum value 
considered. In the present analysis the convergence will be studied with respect to the index 
Kmax, therefore, the number of Laguerre polynomials at each step, rrimax + Ij will be suffi- 
ciently large to guarantee independence from /3 of the physical eigenvalues and eigenvectors. 
We found that rrimax + 1 ~ 20 Laguerre polynomials (with proper values of (3) were sufficient 
for an accuracy of 0.1% in the calculated eigenvalues. 



A. Symmetries of the eigenvectors 

Fixing the total angular momentum and parity J'^ of the state we want to describe, 
the diagonalization of the Hamiltonian produces eigenvectors with well-defined-permutation 
symmetry. Since we are using a central potential, the total angular momentum L and total 
spin S are good quantum numbers. Accordingly, our basis is identified by {L,S,T)J'^, 
where T is the total isospin of the state, and the parity corresponds to consider even or 
odd K values in the expansion. The eigenvalues appear either in singlets, corresponding to 
symmetric or antisymmetric eigenvectors, or in multiplets, corresponding to mixed symmetry 
eigenvectors. The identification of the symmetry of each eigenvector can be done applying 
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to it the Casimir operator 

CiA) = Y,Pihj), (46) 

i<j 

where P(i,j) is the permutation operator of particles Using the results of the preceding 

section, the representation of the Casimir operator in the HH basis resuhs 

C{A) = J2b^''{-1)'^-B^''. (47) 



i<j 



As discussed in Ref. 



23| . this Casimir operator corresponds to the class sum [(2)]^^ of the 
group of permutation of A objects Sa and the corresponding eigenvalues A for the different 
symmetries [A] are given in that reference up to A = 5. The eigenvectors of the Hamilto- 
nian are also eigenvectors of this Casimir operator, therefore the application of the Casimir 
operator to a specific eigenvector \1/^''([A]) results 

C(^)*r([A]) = Av[/f ([A]). (48) 

The different symmetries characterizing the spatial eigenvector are identified by A. The 
physical state of A nucleons is obtained after multiplying \1/^'"([A]) by the proper spin-isospin 
state in order to obtain an antisymmetric state. 



B. A = 3,4: systems 



In Ref. [10| the binding energies -£^0(3) and -E'o(4) corresponding to the ground states 
of the A = 3,4 systems has been studied using the Volkov potential. Here we extend the 
analysis to some more states of the spectrum using, in addition, the s-wave version of the 
potential. In particular the A = 3 system present a very shallow excited state. This is 
a consequence of the very shallow two-nucleon binding energy E2N and the large value of 
the scattering lenth a2N that this potential produces. It is known that when the two-body 
system presents these characteristics, the three-body system could show a certain numbers of 



bound states close to the two-body threshold called Efimov states (see Ref. [2J] and reference 
therein). In the present case, this behavior is a consequence of the parametrization of the 
Volkov potential that has been tuned to approximate the binding energy of the A = 3 
system. In doing that, the binding energy of the A = 2 system results to be much lower 
than the experimental deuteron binding energy. Despite this unrealistic situation, here we 
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are interested in studying the HH expansion for systems with A > 4. The analysis of the 
A = 3,4 systems serves as a basis for estabhshing the different thresholds that appear in 
the description of those systems. 

In Table [T] the A = 3 results for the state (L, 5,T)J'^ = (0, 1/2, 1/2)1/2+ are given using 
the complete potential as well as its s-wave version. The ground state binding energy Eq 
converges at the level of 0.1 keV with K^ax = 40 and, fixing the non-linear parameter 
/3 = 2 fm~^, with rrimax = 24. For the sake of comparison the results of the stochastic 



variational model (SVM) of Ref. 21| and those from Ref. |8|] are given in the table. The 
convergence of the binding energy Ei of the shallow state at the same level of accuracy 
necessitates a much larger basis. The maximum grand angular quantum number has been 
increased up to K^ax = 320 and, with (5 = 1 fm~^, the maximum degree of Laguerre 
polynomials used was rrimax = 32. Above Kmax = 60 only symmetric states with li = I2 = 
have been considered. This very different pattern of convergence in the two binding energies, 



Eq and Ei, has been observed before 25|, |26|]. Moreover, the pattern of convergence of the 



all-waves and s-wave potentials is similar. Since the structure of these states corresponds 
mostly to have the particles in a relative / = state, there is only a small decrease in 
energy for Eq, of about 35 keV, when the s-wave potential is considered. In the excited 
state El this difference is even less, of about 0.5 keV, giving both versions of the potential 
very close values. This is a manifestation of the particular structure of the Efimov state 
in which the third particle orbitates around the / = state of the other two in a very far 
orbit. Increasing the attraction of the two-body potential, the two-body binding energy 
E2N increases faster than the Ei energy and, at so me p oint, the Efimov state starts to be 



above the two-body threshold (see for example Ref. 27|). When realistic forces are used to 
describe the three-nucleon system there is no observation of an excited state, in agreement 
with the experimental situation. However the effective range function presents a pole close 



to the two-body threshold [28|], that can be interpreted as an Efimov-like state embedded in 
the continuum. When the Coulomb interaction is included, the ground state binding energy 
results Eq = 7.7594 MeV (all- waves potential) and Eq = 7.7254 MeV (s-wave potential). In 
both cases the isospin components T = 1/2 and 3/2 are automatically included. With the 
repulsion induced by the Coulomb potential the excited state is not any more bounded. 

The L = state of the A = 4 system is firstly analyzed. The spatially symmetric state 
of four nucleons can be antisymmetrized using the S* = 0, T = spin-isospin functions. 
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In Table ITTl the pattern of convergence, in terms of K^ax, is shown for the first two levels 
of the {L,S,T)J'^ = (0,0,0)0"'" state using both versions of the Volkov potential. The 
ground state binding energy Eq converges at the level of 1 — 2 keV for Kmax = 40 whereas 
the convergence of the excited state binding energy Ei has been estimated at the level of 
50 keV. For both types of potentials the excited state results to be bounded with respect to 
the 3 + 1 threshold. For the sake of comparison the results of the SVM and those from Ref. 8| 
are shown in the table. In order to compare the results to the experimental value of the a- 
particle, i?('^He) = 28.30 MeV, the last four columns of the table show the results including 
the Coulomb potential between the two protons. The obtained values of 29.60 MeV (all- 
waves potential) and 29.43 MeV (s-wave potential) show a pronounced overbinding. This is 
the usual situation when central interactions are used to describe the ^He nucleus and it is 
at variance to the case in which realistic NN forces are used. When the Coulomb potential is 
included the excited state appears slightly above the 3 + 1 threshold. In the case of the all- 
waves potential the lowest threshold, corresponding to a p-^B. configuration, is at 8.465 MeV 
whereas for the s-wave potential it results to be at 8.431 MeV. The n-'^He thresholds are at 
7.759 MeV and 7.725 MeV respectively. Though the convergence for Ei was not completely 
achieved, the description is close to the experimental observation of a 0"'" resonance between 
both thresholds and centered 395 keV above the p-^B. threshold. 

Let us consider the negative parity L = 1 state. The lowest level corresponds to the 
[3 1] irreducible representation and can be antisymmetrized using the S* = 1,T = or 
S = 0,T = 1 spin-isospin functions of four nucleons. Accordingly, using a central potential, 
the J'^ = 0~,1~,2~ states are degenerated. The results are given in Table [XTTl We can 
observe that the all- wave potential produces a bound state at approximate 10.4 MeV far from 
the experimental observation of a 0~ resonance 800 keV above the O"'" resonance. Conversely, 
using the s-wave potential the level results to be unbounded. It appears at approximate 
1.4 MeV above the 3+1 threshold and at approximate 1.3 MeV above the O'^ resonance in 
better agreement with the experimental situation. When the Coulomb potential between 
the two protons is considered the triple degeneracy of the [3 1] representation breaks in 
three different levels, Eq, Ei, E2, showed in the last three columns of the table. The state 
corresponding to the Eo level is formed by an antisymmetric proton pair times a symmetric 
neutron pair and can be completely antisymmetrized with the spin state of the two protons 
having Sp = 1 and the spin state of the two neutrons having Sn = 0, having total spin 
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S = 1. The state corresponding to the Ei level is formed by a symmetric proton pair 
times a symmetric neutron pair and can be completely antisymmetrized with the spin of 
the two protons Sp = and the spin of the two neutrons S'„ = 0, having total spin S* = 0. 
Finally, the state corresponding to the E2 level is formed by a symmetric proton pair times 
an antisymmetric neutron pair and can be completely antisymmetrized with the spin of the 
two protons Sp = and the spin of the two neutrons S'„ = 1, having total spin S = 1. The 
first and third level are mostly T = and can be identified with the {J'",T) = (0~,0) and 
(2~, 0) resonances whereas the Ei level is mostly T = 1 and can be identify with the (1^, 1) 



29|. 



resonance, observed in the low energy spectrum of ^He 

We can conclude that besides its simplicity, the s-wave potential describes the A = 4 
system better than the complete potential and, in some cases, in reasonable agreement with 
the experiment. From a technical point of view we were able to describe the L = ground 
and first excited states and the first level of the L = 1 state using the non-symmetrized HH 
functions. In particular the convergence of the L = first excited state, Ei, presents some 
difficulties since its energy results to be very close to the threshold. When the Coulomb 
interaction is taken into account this level moves to the continuum and it appears as a 
resonance between the two 3 + 1 thresholds, in agreement with the experimental data. The 
s-wave potential describes better also the negative parity resonances; moreover, in order to 
accurately extract their position and width, the present method can be combined with the 



procedure developed for example in Ref. 
in Figs. 



30|. The computed states for A = 3, 4 are collected 



C. A = 5,6 systems 

In the case of systems with A > 4 the spatially-symmetric state cannot be antisym- 
metrized using the corresponding spin-isospin functions. Therefore, it is interesting to study 
the symmetry of the different levels in the A = 5 system when the non-symmetrized basis 
is used. For the positive parity L = state, using the Volkov potential, we found that 
the deepest two levels correspond to a completely symmetric state (the irreducible repre- 
sentation of 5*5 [5]) as expected. As mentioned, they cannot be antisymmetrized using the 
spin-isospin functions of five particles and, therefore, they do not represent physical states 
for five nucleons. The third level belongs to the irreducible representation of 6*5 [4 1]; it 
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can be antisymmetrized using the A = 5 spin-isospin functions having 5* = 1/2, T = 1/2, 
and accordingly it represents the lowest level of the (L, S,T)J'" = (0, 1/2, 1/2)1/2+ state of 
five nucleons. The convergence of these three states in terms of Kmax is given in Table IIVI 
The first two levels, representing bosonic bound states, present a good convergence with K 
(in particular the deepest level). The convergence of the [4 1] state shows that it does not 
describe a bound state, in agreement with the fact that the A = 5 nucleus does not exist. 
In fact its energy results to be above the threshold of 30.42 MeV describing an ^He nucleus 
plus a fifth nucleon far away (here the Coulomb interaction has not been included). For the 
three levels, their stability as a function of the non-linear parameter /3 is shown in Fig. |H 

The negative parity L = 1 state corresponds to the {L,S,T)J'^ = (1,1/2,1/2)1/2" 
and (1,1/2,1/2)3/2" states which are degenerate using the Volkov potential. Its deepest 
level cannot be spatially-symmetric, and in fact it belongs to the [4 1] representation; as 
before, it can be antisymmetrized using the S = 1/2, T = 1/2 spin-isospin functions of five 
nucleons. In Table IVl the convergence of this level is shown in terms of Kmax for the all- waves 
potential as well as its s-wave reduction. From the table we can observe that the all-waves 
version of the potential predicts a very deep bound state, at 43.03 MeV, whereas the s-wave 
reduction does not. Using the s-wave potential, the A = 5 systems results to be unbounded 
in agreement with the experimental observation. From this analysis we can conclude that 
the fact that the A = 5 nucleus does not exist is the result of a delicate balance between 
the Pauli principle, the short range character of the NN interaction and its weakness in 
p-waves. The s-wave potential, used in the present analysis, represents the extreme case in 
which the interaction in p-waves is considered zero. For these two levels, their stability as 
a function of the non-linear parameter /3 is shown in Fig. H] as the solid line (all-waves) and 
long dashed line (s-wave) respectively. In the last three columns of the table, the three levels 
obtained considering the s-wave Volkov potential plus the Coulomb interaction between two 
protons are shown. The inclusion of the Coulomb interaction breaks the degeneracy of the 
quartet- [4 1] state, producing three different states that can be identified by the residual 
5*2 (S> 5*3 symmetry of the two-protons and three- neutrons sub-systems. The lowest two states 
Eq and El belong to the [1^] (g) [3] and [2] ® [3] representations of 5*2 ® S3 and cannot be 
antisymmetrized with respect to the three neutrons. The third level, £'2, is a doublet state 
since corresponds to a mixed symmetry of the three neutrons and it is symmetric in the two 
protons. Its belongs to the [2] ® [2 1] representation and it can be antisymmetrized with the 
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spin state of the two protons having Sp = and the spin state of the three neutrons having 
Sn = 1/2. Physically this state is describing a scattering state between a neutron and an 
a-particle in = l/2~ and 3/2^. In the present study we are hmiting the description 



to bound states, however, using the method described in Ref. [Sjj it would be possible to 
compute phase-shifts using the L = 0, 1 bound-like states. The extension of the method to 
describe scattering states is in progress. 

In the case of the A = 6 system we concentrate the analysis in the (L, S, T)J'^ = (0, 0, 1)0"'" 
and (0, 1,0)1"^" states. Using a central potential, and disregarding the Coulomb interaction, 
these two states are degenerate. Including the Coulomb interaction between two protons, 
the first state has the quantum numbers of ^He. A direct diagonalization of the six body 
Hamiltonian using the non-symmetrized HH basis, with the Volkov potential, produces a 
spectrum in which the first two levels belongs to the [6] irreducible representation of Sq. 
They are completely symmetric and cannot be antisymmetrized using the A = 6 spin-isospin 
functions. The third level belongs to the [5 1] representation and it cannot be antisymetrized 
too. The fourth level belongs to the [4 2] representation, and it is the first one that can be 
symmetrized using the A = 6 spin-isospin functions having S' = 0,T = lorS'=l,T = 0. 
The convergence pattern of these four levels in terms of K^ax are shown in Table IVTl indicated 
by Ei, i = 1, ... ,4. Similar to the A = 5 case, the Volkov potential acting in all waves 
predicts large binding energies. In particular the binding energy of the physical state results 
to be ~ 67 MeV. Using the s-wave potential a much more reasonable value of ~ 34 MeV 
is obtained for this level. The corresponding convergence is shown in the last column of 
Table IVII indicated by E'^. It should be noticed that in the computation of the spectrum 
using the s-wave potential the E^ is not anymore the fourth level. Other levels belonging 
to the [6] and [5 1] representation gain more energy than the [4 2] level, making difficult 
its correct identification. However, it is possible to restrict the search of the eigenvectors 



to those having a particular symmetry using a symmetry-adapted Lanczos method [12| (a 
description of the iterative method used is given in the Appendix). Essentially, starting with 
a vector having the desired symmetry, after each iteration of the matrix- vector product, the 



12|, 



new vector is projected onto the sub-space of the selected symmetry. Following Ref. 
an intermediate purification step is also implemented. This method has the characteristic 
of finding eigenvalues corresponding to eigenvectors of one particular symmetry simplifying 
the search procedure and the identification of the eigenvectors. 
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When the Coulomb interaction between two protons is considered the degeneracy of the 
[4 2] level (of dimension 9) is broken and four different states appear. It is possible to identify 
the physical state looking at the symmetry of the four neutrons. One of the states belongs to 
the [4] representation, two belong to the [3 1] representation and the last one belongs to the 
[2^] representation of S4. This last state is the only one that can be antisymmetrized using 
the spin functions of four neutrons having Sn = 0. Moreover, the proton state is spatially 
symmetric and therefore can be antisymmetrized with the spin function Sp = making a 
total S = state. The convergence of this state is given in the last column of Table I VII It 
should be noticed that this state is embedded in a very dense spectrum. In order to follow 
this state in the projected Lanczos method a projection-purification procedure is performed. 
Essentially the vector, after each matrix vector product, is projected on antisymmetric 
states between particles (3,4) and (5,6). In this way the level belonging to the [2] (g) [2^] 
representation of 5*2 ® 5*4 results to be the lowest one. The results obtained for the different 
levels are collected in Figs. I2]l3l 



V. CONCLUSIONS 

In this work we have developed a technique devoted to describe bound states in an A- 
body system without imposing a particular requirement due to the intrinsic statistic of the 
particles. However, the final aim of the method is to found wave functions that fulfill this 
requirement. 

Starting with the non-symmetrized HH basis set, we have diagonalized the Hamiltonian 
of the A-hodj system using that basis at fixed values of K. We have observed that the 
eigenvectors reflect the symmetries present in the Hamiltonian and, in particular, if the sys- 
tem is composed by identical particles, the eigenvectors belong to the different irreducible 
representations of the permutation group of A objects, Sa- Using a Casimir operator, it 
was possible to identify those eigenvectors having the required symmetry of the system and, 
accordingly, study the convergence (in terms of K) of the corresponding eigenvalues. The 
direct use of the non-symmetrized HH basis has important consequences from a technical 
point of view. The size of the basis is much bigger than the one limited to a subspace having 
a particular symmetry. However, it should be noticed that a system of nucleons includes 
spatial, spin and isospin degrees of freedom, all of them coupled by the NN potential, with 
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the consequence that different spatial symmetries are present in an A-nucleon wave func- 



tion. Although the construction o 



HH basis elements having different spatial symmetries 



is possible (see for example Ref. js]), the necessity of including the different symmetries in 
the description enlarges the dimension of the basis and makes it comparable to the case in 
which the non-symmetrized basis is used. This is particularly important when one wants to 
consider the description of the small components of the wave function induced by symmetry 
breaking terms in the potential, as for example high isospin components. 

The method here presented is based in a particular implementation of the potential 
energy matrix constructed as a sum of products of sparse matrices. This allows to efficiently 
use iterative algorithms in which the matrix- vector product is a key element. However the 
iterative methods are well suited to calculate the deepest levels of the spectrum. In our 
formulation, due to the presence of different symmetries, the physical states could appear 
very high in the spectrum or in a zone with a high density of levels. In this case we found 



very convenient to use the symmetry-adapted Lanczos method [12|. Using the particular 
form of the permutation operator P{i,j) in terms of the sparse matrices (see Eq. (l47|) ). it was 
possible to project the vector in the iterative procedure to be antisymmetric in selected pairs 
of particles. In this way the desired symmetry becomes the lowest state of the spectrum. 
Though this mechanism is not as fast as searching for the true lowest state of the complete 
spectrum, it is much faster than searching for certain numbers of levels in high position of 
the spectrum. 

We should also stress that the sparse matrices ^^^A defined in Eq. (l30|) have the property 
of being constructed as products of the angular T-coefficients, the tree T-coefficients and the 
Raynal-Revai coefficients. Whereas the latters couples quantum number belonging to the [K] 
and [K'] sets, the T- and T-coefficients perform a recoupling of quantum numbers inside [K] 
or [K']. Moreover the Raynal-Revai coefficients T^f^i^ini couple quantum numbers belonging 
to three particles (see Eq.( p8|) ). As the number of particles A increases, more values of the 
quantum numbers L are accessible and this makes the size of the '^'^'^A matrices to increase, 
though slowly with the number of particles A. Furthermore, going from a system with A 
particles to A + 1, the number of potential terms increase by A and the number of factors in 
the matrix Bf^^' of Eq.(l33]) increases at maximum of the same quantity. The computational 
effort increases roughly linear with A and this fact makes feasible the application of the 
method for increasing values of A as has been demonstrated in the present work. Our 
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expectation is that the present technique could be extended to treat systems up to A = 8. 
The calculations presented here have been obtained using a sequential code. We expect that 
an opportune parallelization of the code (which is under study) will increase the potentiality 
of the method. 



We have limited the analysis to consider a central potential, the Volkov potential, used 
several times in the literature. Though the use of a central potential leads to an unrealistic 
description of the light nuclei structures, the study has served to analyze the characteristic 
of the method: the capability of the diagonalization procedure to construct the proper 
symmetry of the state and the particular structure, in terms of products of sparse matrices, 
of the Hamiltonian matrix. The success of this study makes feasible the extension of the 
method to treat interactions depending on spin and isospin degrees of freedom as the realistic 



NN potentials. A preliminary analysis in this direction has been done [32|. To this respect 
it is important to notice that the information of the potential is given in the matrix Vu 
defined in Eq. fl36l) . Once this matrix is given, the method remains the same, with the 
basis enlarged to include spin and isospin degrees of freedom if required. Using the Volkov 
potential we have shown that it was possible to identify all the physical states and the 
corresponding thresholds in order to interpret the level as bounded or belonging to the 
continuum. Furthermore, the results obtained using the Volkov potential up to A = 6 
compare well with other techniques. A few characteristic of the A = 3 — 6 systems using 
the Volkov potential are the following. Due to its particular parametrization a shallow state 
appears in the A = 3,4 systems when the Coulomb interaction is not considered. In the 
A = 3 this state has the characteristic of an Efimov state. When the Coulomb interaction 
is considered these states move to the continuum. The Volkov potential acting in all waves 
produces large binding energies as A increases. Accordingly we have included in the analysis 
the s-wave version of the potential. In agreement with the experimental observations, this 
version predicts in reasonable positions the A = 40~^ and 0~ resonances and no bound states 
in the A = 5 system. It also predicts reasonable binding energies in the A = 6 system. The 
extension of the method to consider realistic potentials is in progress. 
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VI. APPENDIX 



The diagonalization of the Hamiltonian is obtained by means of an iterative algorithm 
which requires only the action of the Hamiltonian matrix on a given vector. We used 



the Lanczos algorithm in the version invented by CuUum and Willoughby [33| which is 
particularly sparing with memory use. In principle, the iterative procedure should preserve 
the permutation symmetry of the input vector, as the Hamiltonian commutes with the group 
elements. However, the round-off errors generate components also in the other irreducible 
representations. To circumvent this problem, we have used a symmetry adapted Lanczos 



(SAL) developed in Ref. 12| in which a projection operator is applied after each iterative 



step. Starting from a random initial vector, in the usual Lancsoz recurrence formula 

Pi+iVi+i = Hvi - aiVi - PiVi^i , (49) 

the product Hvi is replaced by P^^^Hvi, where P'"^' is a projector on a sub-space with a 
non-zero intersection with the irreducible representation [A], and zero intersection with the 
irreducible representations of the lower-eigenvector symmetries. A purification step is also 
performed in which the product (3i+iVi^i is replaced by Pl'^l/^j+i'Uj+i. 

As an example in the L = sector of the A = 6 system, we are interested in states 
belonging to the irreducible representation [4 2]. In order to eliminate lower states belonging 
to the irreducible representations [6] and [5 1], we have used the projector 

P[^]=A12-A34, (50) 

given as the product of the antisymmetrization operator A12 with respect particles (1,2), and 
the antisymmetrization operator ^34 with respect particles (3,4). The two antisymmetriza- 
tion operators have the following expression in terms of the Ai matrices (the superscript 
L, M = 0, is understood) 

^12 = ^(1-^5), (51) 



and 

A34 = AaAsAa ^""^ A4A3A4 . (52) 



1-A, 



25 



[1] S.C. Pieper, K. Varga, and R.B. Wiringa, Phys. Rev. C66, 044310 (2002). 

[2] P. Navratil, V.G. Gueorguiev, J. P. Vary, W.E. Ormand, and A. Nogga, Phys. Rev. Lett. 99 
042501 (2007). 

[3] A. Kievsky et al., Phys. Rev. CSS, 3085 (1998). 

[4] R. Lazauskas et al., Phys. Rev. C71, 034004 (2005). 

[5] H. Kamada et al., Phys. Rev. C64, 044001 (2001). 

[6] A. Novoselsky and J. Katriel, Phys. Rev A 49, 833 (1994). 

[7] A. Novoselsky and N. Barnea, Phys. Rev A 51, 2777 (1995). 

[8] N. Barnea, Phys. Rev. A 59, 1135 (1999). 

[9] N. K. Timofeyuk, Phys. Rev. C 78, 054314 (2008). 

[10] M. Gattobigio, A. Kievsky, M. Viviani, and P. Barletta, Phys. Rev A 79, 032513 (2009). 

[11] M. Viviani, A. Kievsky, and S. Rosati, Phys. Rev. C 71, 024006 (2005). 

[12] Xiao-Gang Wang and Tucker Carrington, Jr, J. Chem. Phys. 114, 1473 (2001) 

[13] N. Ya. Vilenkin, G. I. Kuznetsov, and Ya. A. Smorodinskii, Sov. J. Nucl. Phys. 2, 645 (1966) 

[14] M.S. Kil'dyushov, Yad. Fiz. 15, 197 (1972) [Sov. J. Nucl. Phys. 15, 113 (1972)] 

[15] M.S. Kil'dyushov, Yad. Fiz. 16, 217 (1972) [Sov. J. Nucl. Phys. 16, 117 (1973)] 

[16] N. Barnea and A. Novoselsky, Ann. of Phys.256, 192 (1997). 

[17] J. Raynal and J. Revai, Nuovo Cimento A 68, 612 (1970). 

[18] R. Krivec and V.B. Mandelzweig, Phys. Rev. A 42, 3779 (1990). 

[19] M. Viviani, Few-Body Syst. 25, 177 (1998). 

[20] V.D. Efros, Few-Body Syst. 19, 167 (1995). 

[21] K. Varga and Y. Suzuki, Phys. Rev. C 52, 2885 (1995). 

[22] N. K. Timofeyuk, Phys. Rev. C 65, 064306 (2002). 

[23] A. Novoselsky and J. Katriel, Phys. Rev C 51, 412 (1995). 

[24] E. Braaten and H.-W. Hammer, Phys. Rept. 428, 259 (2006) 

[25] P. Barletta and A. Kievsky, Few-Body Syst.45, 25 (2009) 

[26] P. Barletta, C. Romero-Redondo, A. Kievsky, M. Viviani, and E. Garrido, Phys. Rev. Lett 

103, 090402 (2009) 

[27] P. Barletta and A. Kievsky, Phys. Rev. A 64, 042514 (2001) 

26 



[28] A. Kievsky et al., Phys. Lett. B 406, 292 (1997) 

[29] D.R. Tilley, H.R. Weller and G.M. Hale, Nuc. Phys. A 541, 1 (1992) 

[30] H. Witala and W. Glockle, Phys. Rev. C 60, 024002 (1999) 

[31] A. Kievsky, M. Viviani, P. Barletta, C. Romero- Redondo, and E. Garrido, Phys. Rev. C 81, 
034002 (2010) 

[32] M. Gattobigio, A. Kievsky, M. Viviani, and P. Barletta, Few-Body Syst. 45, 127 (2009) 

[33] J. K. Cullum and R. A. Willoughby, J. Comput. Phys. 44, 329 (1981) 



27 



TABLE I. yl = 3 results for (L, S, T)J'' = (0, i, \)\ state usine the all-waves and s-wave Volkov 
potential as a function of the maximum grand angular quantum number Kmax- The ground state 
Eq as well as the excited state Ei are shown. 





all- waves 


s- 


wave 


K 


Eq (MeV) El (MeV) 


Eo (MeV) El (MeV) 


20 


8.4623 0.3627 


8.4283 


0.3618 


40 


8.4649 0.5181 


8.4309 


0.5174 


60 


8.4649 0.5595 


8.4309 


0.5589 


80 


8.4649 0.5773 


8.4309 


0.5768 


100 


0.5866 




0.5861 


120 


0.5918 




0.5913 


140 


0.5947 




0.5943 


160 


0.5965 




0.5960 


180 


0.5976 




0.5971 


200 


0.5982 




0.5978 


240 


0.5989 




0.5985 


280 


0.5992 




0.5988 


320 


0.5993 




0.5989 


SVM [21] 


8.46 






Ref. [8] 


8.462 0.2599 
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TABLE II. Binding energies for the A = 4 ground state £"0 and the first excited state Ei of the 
(L, S, T)J'" = (0, 0, 0)0''" state using the ah-waves and s-wave Volkov potentials as a function of the 
maximum grand angular quantum number K^ax ■ In the last four columns the Coulomb interaction 



has been considered. For the sake of comparison the results of Refs. 



2l| are shown. 



all-waves s-wave 
Eq (MeV) El (MeV) Eq (MeV) Ei (MeV) 






28.580 


3.238 


28.580 


3.238 


27.748 


2.787 


27.748 


2.787 


10 


30.278 


7.509 


30.116 


7.445 


29.456 


7.039 


29.292 


6.976 


20 


30.416 


8.223 


30.250 


8.164 


29.596 


7.778 


29.429 


7.720 


30 


30.418 


8.463 


30.252 


8.403 


29.599 


8.035 


29.431 


7.976 


40 


30.418 


8.562 


30.252 


8.501 


29.600 


8.144 


29.432 


8.085 


SVM [21] 


30.42 
















Ref. [8] 


30.406 


8.036 















all-waves s-wave 
£^0 (MeV) El (MeV) Eq (MeV) Ei (MeV) 



TABLE III. The binding energy of the A = 4 lowest level having L = 1, using the all- waves and 
s-wave Volkov potential, as a function of the maximum grand angular quantum number Kmax- In 
the case of s-wave potential, when Coulomb interaction between particles (1,2) is considered, the 
level splits in three sub-levels, whose energies Eq, Ei and E2 are shown in the last three columns. 



Kfnax 


all- waves 


s-wave 


Eq (MeV) El (MeV) E2 (MeV) 


1 


7.965 


0.387 








3 


8.411 


1.975 


1.639 


1.440 


1.374 


11 


10.121 


5.567 


5.314 


5.091 


4.899 


21 


10.373 


6.642 


6.456 


6.276 


5.955 


31 


10.406 


7.113 


6.965 


6.850 


6.417 
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TABLE IV. A = 5 binding energies of the first tliree levels of the L = state, belonging to the 
indicated irreducible representation [A], as a function of K^ax- The size of the HH basis Nhh is 
also indicated. 



K 


Nhh 


Eq (MeV) El (MeV) E2 (MeV) 






[5] 


[5] 


[4 1] 





1 


64.864 


24.472 


- 


2 


10 


64.864 


24.472 


20.160 


4 


55 


65.958 


28.411 


22.043 


6 


220 


66.893 


29.517 


24.415 


8 


714 


67.713 


30.228 


25.568 


10 


1992 


68.008 


30.587 


26.459 


12 


4950 


68.177 


30.927 


27.043 


14 


11220 


68.239 


31.152 


27.515 


16 


23595 


68.264 


31.357 


27.862 


18 


46618 


68.274 


31.509 


28.143 


20 


87373 


68.278 


31.628 


28.371 


22 


156520 


68.279 


31.715 


28.560 


24 


269620 


68.280 


31.779 


28.719 
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TABLE V. ^ = 5 binding energies of the deepest L = 1 state, as a function of K^axi using the 
all-waves and s-wave Volkov potential. In the last three columns the Coulomb potential has been 
summed to the s-wave Volkov potential. The size of the HH basis Nhh is also indicated. 



Kmax 


Nhh 


all-waves 


s-wave 


Eq 


El 


E2 


1 


4 


39.635 


21 


874 


21.370 


21.119 


- 


3 


40 


40.001 


24 


317 


23.854 


23.604 


23.524 


5 


220 


41.022 


26 


053 


25.618 


25.367 


25.251 


7 


876 


41.785 


26 


923 


26.505 


26.258 


26.116 


9 


2820 


42.384 


27 


546 


27.140 


26.896 


26.736 


11 


7788 


42.682 


27 


971 


27.574 


27.333 


27.160 


13 


19140 


42.868 


28 


297 


27.908 


27.669 


27.485 


15 


42900 


42.952 


28 


521 


28.140 


27.903 


27.710 


17 


89232 


42.996 


28 


693 


28.320 


28.084 


27.882 


19 


174460 


43.017 


28 


823 


28.457 


28.223 


28.011 


21 


323752 


43.027 


28 


924 


28.562 


28.331 


28.110 


23 


574600 


43.032 


29 


005 


28.647 


28.417 


28.189 


SVM 
HH [8j 




43.00 
42.383 
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TABLE VI. A = 6 binding energies of the first four levels of the L = state, using the Volkov 
potential, belonging to the indicated irreducible representation [A], as a function of Kmax- 
indicates the binding energy of the lowest [4 2] state using the s-wave Volkov potential, and E^'^ 
is the binding energy of the [2] [2^] state, once the Coulomb interaction has been included. The 
size of the HH basis Nhh is also indicated. 

Kmax Nhh Eq (MeV) Ei (MeV) E2 (MeV) E^^ (MeV) ^| (MeV) Ef (MeV) 
[6] [6] [5 1] [4 2] [4 2] [2]® [2^] 






1 


117, 


.205 


64, 


.701 












2 


15 


117, 


.205 


64, 


.701 


62.513 


61.142 


24. 


,793 


24.064 


4 


120 


118, 


.861 


69, 


.450 


64.277 


62.015 


28. 


,791 


28.016 


6 


680 


120, 


.345 


70, 


.544 


66.268 


63.377 


30. 


,723 


29.935 


8 


3045 


121, 


.738 


71, 


.443 


67.280 


64.437 


31. 


,645 


30.851 


10 


11427 


122, 


.317 


71, 


,923 


68.371 


65.354 


32. 


,244 


31.446 


12 


37310 


122, 


.597 


72, 


,477 


69.029 


65.886 


32. 


,708 


31.908 


14 


108810 


122, 


.711 


72, 


.822 


69.531 


66.201 


33. 


,075 


32.275 


16 


288990 


122, 


.752 


73, 


,101 


69.842 


66.360 


33. 


,358 


32.558 


18 


709410 


122, 


.768 


73, 


,284 


70.051 


66.437 


33. 


,561 


32.762 


20 


1628328 


122, 


.774 


73, 


.407 


70.189 


66.474 


33. 


,710 


32.912 


22 


3527160 


122, 


.776 


73, 


.485 


70.283 


66.491 


33. 


,814 


33.017 


SVM 














66.25 
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[2] 0+ 0.599 [3] 0+ 
A = 2 



8-465 [3] 0+ 8.562 [4] 0+ 

^ = ^ 1 0.406 [31] 1- 

28.72 [41] 0+ 

30.418 [4] 0+ 

31.72 [5] 0+ 



A = 4 



43.03 [41] 1- 

66.49 [42] 0+ 
68.28 [5] 0+ — — 

70.28 [51] 0+ 



73.49 [6] 0+ 
122.78 [6] 0+ 



FIG. 2. Calculated levels for A = 2 — 6 using the all-waves Volkov potential. 
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0.546 0+ 


0.599 0+ 








9tt 










6.417 2-,0 
6.850 1-,1 




(MeV) 




8.431 0+ 


6.965 0-,0 

7.725 0+ 

8.085 0+,0 
^He — 




Energy 






2Q 4S 0+ 


33.02 0+ 










^He 



FIG. 3. Calculated levels for A = 2,3,4 and 6, using the s-wave Volkov potential with the inclusion 
of Coulomb interaction for He isotopes. In this case the A = 5 system results unbounded. 



35 



> 

CD 



-25 
-30 
-35 
-40 
-45 



^ -50 
W -55 



-60 
-65 
-70 







L = 1 all- waves 

L = 1 s-wave 

L = Eq 

L = El 

L = E2 

_i I I I 

12 3 4 

^ (fm-i) 



FIG. 4. The A = 5, L = levels, given in Table ITVl denoted as Eq, Ei, E2 and the L = 1 levels given 
in Table |Vl denoted as all-waves and s-wave, are shown as functions of the non-linear parameter 
(3, at Kmax = 16 (L = levels) and Kmax = 17 (L = 1 levels), respectively. 
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